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Abstract 

Taking precautions before or during the start of a virus outbreak can heavily reduce the number of 
infected. The question which individuals should be immunized in order to mitigate the impact of 
the virus on the rest of population has received quite some attention in the literature. The dynamics 
of the of a virus spread through a population is often represented as information spread over a 
complex network. The strategies commonly proposed to determine which nodes are to be selected 
for immunization often involve only one centrality measure at a time, while often the topology of 
the network seems to suggest that a single metric is insufficient to capture the influence of a node 
entirely. 

In this work we present a generic method based on a genetic algorithm (GA) which does not rely 
explicitly on any centrality measures during its search but only exploits this type of information to 
narrow the search space. The fitness of an individual is defined as the estimated expected number 
of infections of a virus following SIR dynamics. The proposed method is evaluated on two contact 
networks: the Goodreau's Faux Mesa high school and the US air transportation network. The GA 
method manages to outperform the most common strategies based on a single metric for the air 
transportation network and its performance is comparable with the best performing strategy for the 
high school network. 

Keywords: Complex networks; Genetic algorithm; Node centrality; SIR model; Damage control 

1 Introduction 

Many viruses spread through human population by means of personal contact between infectious indi- 
viduals (those who carry the virus) and susceptibles (those who are not ill at the moment but can catch 
the disease) HH3. A concept that proved to be very valuable in order to gain a better understanding 
of the virus spread process is the complex network [3]. Nodes within this graph structure represent 
individuals and are associated with a certain state (e.g., susceptible or infectious). Edges between these 
nodes account for social interactions. The percolation of the virus through the population (network) is 
then predicated by a fixed set of rules [2j|4j|5l|6l. 

This approach has not only been employed in order to better understand the dynamics of such a 
disease spread, but also gives rise to a field of research with a more proactive attitude: which individuals 
in the population should be immunized to limit the spread of the virus most effectively? Or, to frame it 
a bit differently, which nodes in the network should be protected in order to limit the damage done as 
much as possible? 

The concept commonly used in order to find those nodes that are in need of protection (or removed 
in some cases) is node centralit^\ Q |U [TO] \TT\ Q21. Node centralities express to what extent the 

* Corresponding author. Tel.: +7-965-073-0861. E-mail address: kashirin.victor@gmail.com. 
'Measures for edge centrality exist as well, see (3). 
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Fig. 1: A rather simple hypothetical network. The values in each node represent the node's degree and 
betweenness centrality, respectively. 



node facilitates the spread of information over the network (note that information spread is, to some 
extent, alike to a virus spread lfT3l ). In this article we consider three common traditional variant^] used 
in the literature, each of which formalizes the concept of centrality in a (slightly) different manner: 1) 
degree centrality which is equal to the node's degree, 2) betweenness centrality which focusses rather 
on to what extent the node could influence the communication between other nodes in the network, and 
3) eigenvector centrality which expresses the relative importance of a node in terms of importance of 
its direct neighbors (a having a few connections with important nodes is more important than the one 
having the same number of connections with less important ones). The definitions of these centrality 



measures are given in section 2.3 The usual approach to locate those nodes that influence the spread of 
information the most is to determine for all nodes in the network their score on one centrality measure 
only. A fraction of nodes that scored the highest are then proposed to be immunized. This approach 
has its drawbacks. Take, for example, the network in Fig. [T] It consists of two clusters connected by 
one single node. If one would only use a ranking based on degree centralities (as it rather common in 
the literature), one would start by immunizing several nodes in the clusters. The node in the center of 
the graph is left untouched, while immunizing it at an early phase, would prohibit the virus to spread 
from one cluster to the other. The importance of this node would be noted only when one would take 
betweenness centralities into account as well. Of course, this example is rather simple and one would be 
able to come up with such an analysis by simply examining the graph visually. When the network gets 



large, however, (e.g., 500 nodes as in one of the networks examined later - see section 2.2) it is rather 
hard to make these kinds of assertions and a more structural approach is required. 

In this paper we propose a genetic algorithm [ 14] (GA) for finding those individuals that need to be 
protected in order to limit the spread of the virus through the rest of the network as much as possible. 
This heuristic search approach differs from the methods proposed in the literature, since the GA is 
allowed to search through a wide range of subsets of nodes and does not rely explicitly on any kind of 
centrality metrics during its search; centrality measures are only used to narrow the search space. Nodes 
that score low on all measures of centrality most likely do not play an important role in the spread of the 
virus and can, thus, be neglected. 

In order to emulate the spread of the virus over the network we define an adapted version of one of 
the most common models in epidemiology: the SIR model fl~l|2l[3j[l3]]. The model was first introduced 
in the twenties by Lowell Reed and Wade Hampton Frost. They proposed to divide the entire population 
into the following three distinct classes: 

• The susceptibles, S: the group of individuals who have not been infected but can catch the disease. 

• The infectives, I: the group consisting of all individuals that are currently infected by the disease 
and could infect others from the susceptible class. 

• The removed, R (sometimes referred to in the literature as recovered): those who had the virus but 
either recovered and gained immunity or died. 



-More complicated centrality measures have been proposed as well, see, for example, 191 1111112 
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Note that an individual in the model can go through three strictly sequential phases: susceptible indi- 
viduals can get infected by t he vi rus and either recover or die; removed individuals never loose their 
gained immunit)^] In section 2.4 we discuss our adaptation of this model for simulating the spread on 
networks. 

The structure of the paper is organized as follows. In the following section we discuss the methods 



used. First we define the optimization task at hand formally in section [271] Section 2.2 introduces the 
two networks that are used for validation of the methods. We proceed with defining the three different 
centrality measures discussed earlier in the text. In section |2.4| we formally define an adaption of the 



traditional SIR compartmental model for network structures. Section|23]contains a detailed description 
of the genetic algorithm for finding the optimal set of nodes that need to be protected using the various 
centrality rankings. The results are presented in section [3j We finish this paper with our conclusions, 
some discussion and a few pointers for future research. 



2 Methods 



2. 1 The optimization task 

The optimization task presented in the introduction can be expressed a bit more formally in the following 
manner. Suppose we are dealing with a network G = (V,E) where V denotes the set of vertices in the 
network and E is the set of (bidirectional) edges between pairs of nodes in V. The size of the network 
is given by N = \V\. The number of nodes that can be protected is limited due to time and resource 
constraints to a total of k nodes. We are, thus, interested in finding a ^-subset of nodes I QV that, when 
protected, will most effectively limit the spread of the virus through the network G. This task is by no 
means easy. The number of possible ^-subsets of nodes is [A which, for most real-world networks, is 
simply too large to brute-force. In addition, one is generally unaware where the virus starts which needs 
to be accounted for when assessing the successfulness of a solution. 

In the following sections we describe a method to find a (near-)optimal ^-subset of nodes. 

2.2 The data sets used for validation 

In this paper we consider two contact networks to validate our approach, one social and one transporta- 
tion network: 

1. The Goodreau's Faux Mesa high school network which is the result of a simulation of a in-school 
friendship network based on data from a high school in the rural western United States [15]. The 
network consists of 147 nodes and 202 undirected edges (representing mutual friendship). See 
Fig. [4^ for a visual representation of the graph. 

2. The US air transportation network [ 16], which consists of 500 nodes (US airports with the highest 
amount of traffic) and 2980 undirected edges (representing air travel connections). See Fig. [4}) for 
a visual representation of the graph. 

Both data sets are freely available and can be found, respectively, at the website of the CASOS group|^] 
and the Cx-Nets websitd^ 

Please, note that the models, methods and measures discussed in the following sections can be 
applied (or extended) to any network structure. The networks chosen here are only used as examples. 

3 Several extensions of the models exist, e.g., the SIRS model where recovered individuals can return to the susceptible 
state. See the paper by H.W. Hethcote (2000) for a splendid overview of these so-called compartmental models. 
4 http://www.casos.cs. cmu.edu/computationaLtools/datasets/external/Goodreau/indexl l.php 
5 http://sites. google.com/site/cxnets/usairtransportationnetwork 
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2.3 Node centrality 

Measures for centrality try to formalize the (relative) importance of a node in the network. In this paper 
we consider the following usual measures: 1) degree, 2) betweenness and 3) eigenvector centrality. 

Degree centrality, Co, formalizes the importance of a node by setting it equal to the number of 
neighbors the node has; the more connections, the more important the node: 

C D (v)=deg(v). (1) 

Betweenness centrality, Cg, tries to capture a different kind of 'importance'. The idea behind this for- 
malism is that the importance of a node depends on to what extent it can influence the communication 
between the other nodes in the network Q. The measure is computed by determining the fraction of the 
shortest paths between each pair of nodes in the network that contains the node v, i.e.: 

C(v)= £ °^>, (2, 

where G st is the number of the shortest paths between the nodes s and t and a s , t (v) returns the number 
of shortest paths that passes through the vertex of interest, v. As one can imagine, this measure for 
centrality is computationally rather hard to determine. In 2001 Ulrik Brandes presented a fast version 
of the algorithm; computation time is upper bounded by We will use this algorithm for 

determining the betweenness centrality throughout this paper. 

Eigenvector centrality, Ce, distinguishes itself from the other measures discussed here, since it takes 
the importance of its neighbors explicitly into account: 

Cb(v) = T V a vw C E (w), (3) 

where a VM is the (v,w)-entry of the adjacency matrix of the network in question. Rewriting this expres- 
sion to matrix form yields 

Ax = Ax, (4) 

where the z'-fh entry of x is equal to Cg(i). In other words, the eigenvector centrality of a node v is the 
v-th entry of the right eigenvector associated with the first eigenvalue A of the adjacency matrix A of the 
network. 

2.4 An adaption of the SIR model for networks 

The SIR model is one of the most often used models for simulating the spread of a virus through a 
population. The original model in terms of a set of ordinary differential equations makes the assump- 
tions that 1) the population size, N, is large and 2) the population is 'fully mixed', meaning that every 
individual has the same number of (randomly picked) connections and have (approximately) the same 
number of contacts at the same time 01 . Of course, these assumptions are rather unsatisfactory and 
can be overcome by extending the SIR model to networks [2j|3l. 

Suppose we have the following network G = (V,E) where the set V denotes the vertices of the 
network and E is the set of bidirectional edges between pairs of vertices in V. Each node v in the set 
of vertices V is associated with either the state S (susceptible), / (infectious) or R (removed) at each 
moment in time, i.e., St(v,f) G {S, I,R}. The initial state of the network is given by {5'o,/o ; ^o}> i-£-> the 
sets of nodes who's state at the start of the simulation, t = 0, are either S,I or R. The dynamics of the 
virus are emulated by updating the state of all the nodes in the network simultaneously while stepping 
forward in time. When the state of node v G V at moment t is susceptible, i.e., St(v,f) = S, the state of 
this node at the next time step is given by 

S«(,, + l):=( ; wiftprcbabilityl-a-^KOI, & 
S otherwise, 
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where j3 is the chance of the susceptible node to be infected by a single infectious neighbor. The function 
.JV V {1) returns the set of infectious neighbors in the direct vicinity of node v. The operator |.| returns the 
cardinality of a set. Note that the chance of not becoming infected decreases exponentially with the 
number of infectious neighbors. 

When node v is at time t infective, its state becomes 

jR with probability 7, 
St(v,f + 1) := < (6) 
I / otherwise, 

where 7, thus, denotes the chance of recovering (or dying) from the virus during one time step. Since we 
assumed that once a node reached the removed state R, it is either indefinitely immune for the virus or 
dead, the state of that node will not change till the end of the simulation. The simulation ends when there 
are no more infective nodes in the networks; the spread of the virus grinded to a halt and the number of 
susceptible and removed nodes in the network will not change. The total number of casualties is equal 
to the number of nodes in state R. 

2.5 Genetic algorithm 

In order to find a (near-optimal) set of nodes that need to be protected from the virus we employ a genetic 
algorithm [14]. First, we compute the degree, betweenness and eigenvector centrality for each node in 



the network (see section 2.3 ). By ranking the nodes according to their centrality scores, we obtain the 



following three ascending rankings: 

^ = {vfi),vf 2) ,...,vf w) }, R B = {v B {l) ,vf 2) ,...,vf N) } and R E = {vf Y) ,vf 2) , . . . ,vf N) }, (7) 

where Rp, Rb and Re sue the rankings based on, respectively, the degree ([T]), betweenness Q and 
eigenvector ([3]) centrality. We limit the search space of the GA by taking into account that nodes with 
low centrality scores on all three centrality measure do most likely not play an important role in spreading 
the virus over the network. We, thus, keep only the first I <N nodes in each ranking, i.e., 

^{vf^vg),...,^)}, R' B = {v B {l) ,vf 2) ,...,v B {l) } and R' E = {vf^vf^, . . . ,vf }. (8) 

The genetic algorithm will only search for an optimal ^-subset of nodes in the union of these reduced 
rankings: 

R'=R' D UR' B UR' E . (9) 
Each individual in the GA is represented by a chromosome with k genes: 

/ = {vi,V2,...,vjt} where V; £ R' . (10) 

The genes, thus, represent the nodes that need to be protected from the virus where we only take into 
account those nodes that score relatively high on at least one of the centrality measures. 

In order to determine the fitness of an individual, firstly we take the original network G = (V,E) 
and remove the nodes present in the chromosome I. (Removing nodes corresponds here to immunizing 
the nodes from the virus). Secondly, in order to account for the fact that one is normally unaware 
where the virus starts, we select one node at random in the new graph to be infectious while keeping 



the other nodes susceptible. Thirdly, we apply the SIR model as described in section 2.4 and store the 
number of removed individuals at the end of the simulation which we denote with iV casua iti es . (Recall 
that the simulation ends when there are no more infectious nodes in the network). Since the selection 
of infectious individuals is a stochastic process we repeat the last two steps m times. The fitness of an 
individual in the GA population is then defined as the expected number of casualties: 

j jit 

fitness(7) = E [tf casualties ] « " I <Laities dD 
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Table 1 : Parameter settings. 


Symbol 


Description 




Value range 


n 
P 

7 


ihe chance of becoming infected by one infectious neighbor during one time step, see eq. 151 
The chance of recovering during one time step, see eq. J6[ 




.3 

.3 


k 
1 

m 


The number of nodes that are to be protected from the virus 

The number of nodes with the highest scores on a centrality measure that should be kept, see section 2.5 
The number of SIR simulation runs for the deteiniining the fitness of an individual in the GA, see eq. |T 


] 


{10,50} 
{100,200} 
100 




The size of the GA population 

The number of generations simulated by the GA 

The size of the tournament used during tournament selection 

The mutation rate for a gene 




100 
100 
4 

l/k 



where A^Lames * s tne number of recovered nodes at the end of the i-th simulation of the SIR model 
applied to the updated graph G. 

In summary, the objective of the GA is find that subset /* that minimizes the expected number of 
casualtie^J 

I* = argmin fitness (/). (12) 

ICR',\I\=k 

The genetic algorithm simulates 100 generations each with a total of 100 individuals. Selection 
for mating is performed by applying tournament selection where the tournament size is set to 4. As 
crossover operator we use an adaptation of uniform crossover. Suppose we selected two parents to mate, 
e.g., Pi = {1,2,3,4} and P2 = {3,4,5,6}. We then concatenate the chromosomes of these parents and 
sort the resulting array. In our example, we find {1,2,3,3,4,4,5,6}. The chromosome of the first child, 
C\, consists of the odd entries of the array; the chromosome of the second child, C2, consists of the even 
entries, e.g., C\ = {1,3,4,5} and C2 = {2,3,4,6}. This approach guarantees the absence of duplicates 
in the new chromosomes. The mutation rate per gene is set to l/k, i.e., each chromosome undergoes 
on average one mutation in one gene every generation. The mutation of a gene entails that its value is 
randomly set to a node in the set R' that is not in the individual's chromosome already. 

The initial population is randomly generated except for three individuals: their chromosomes are 
set to the first k nodes in, respectively, the degree, betweenness and eigenvector centrality rankings as 
given in eq. ([8]). In addition, we keep for every generation the 10 best performing individuals from the 
previous generation in order not to loose good solutions that were found earlier. 



3 Results 

Table [T] provides an overview of the parameters used in the SIR model and the GA introduced earlier. 
In addition, it presents the parameter settings that were used for producing the results presented in this 
section. 

Fig. [2] depicts the correlations between the degree, betweenness and eigenvector centralities for the 
nodes in both the networks. The reported R 2, s are determined by applying linear regression. Note that 
in both cases degree centrality seems to correlate with its betweenness and eigenvector counterparts 
while betweenness and eigenvector centralities do not seem to correspond. It is interesting to see that 
degree and eigenvector centralities seem to correlate more for the air transportation network (e) than for 
the high school network (b). This might be due to the higher number of edges in the air transportation 
network (2980 against 500 nodes) in contrast to the high school network (202 edges against 147 nodes). 
The main point of these figures is to show that there are discrepancies between the various centralities 



measures, which the genetic algorithm from section 2.5 is set out to exploit in order to find a more 
optimal set of nodes. 



Although it is more common to maximize a fitness function rather than minimizing it, we felt minimization would be more 
appropriate, since we are dealing with the expected number of casualties. 
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Fig. 2: The correlations between the degree, betweenness and eigenvector centralities in both networks. 
The scatter plots (a)-(c) represent the data from the Goodreau's Faux Mesa high school network, and the 
scatter plots (d)-(f) depict the data from the air transportation network. 



So, to what extent does combining centralities measures help in finding a more optimal set of nodes 
that need to be protected from the virus? In order to answer this question, we compare four different 
strategies: 

1. Provide protection for the nodes with the highest degree centrality, see eq. ([TJ; 

2. Provide protection for the nodes with the highest betweenness centrality, see eq. ([2]); 

3. Provide protection for the nodes with the highest eigenvector centrality, see eq. Q; 
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GA degree betweenness eigenvector no protection 



degree betweenness eigenvector no protection 



Fig. 3: The boxplots for the number of infected nodes found when various strategies are applied. GA 
stands for the genetic algorithm, 'degree', 'betweenness' and 'eigenvector' stands for the cases when 
only the nodes with the highest score on that particular centrality measure where protected from the 
virus, 'no protection' refers to the case when no preparations were undertaken. The whiskers represent 
the lowest datum still within 1.5 IQR of the lower quartile and the highest datum still within 1.5 IQR of 
the upper quartile. (IQR stands for interquartile range). Outliers are denoted by a cross. 
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(a) (b) 

Fig. 4: The topology of the Goodreau's Faux Mesa high school network (a) and the US air transportation 
network (b). The nodes depicted with A are the nodes that were suggested by both the genetic algorithm 
and the strategy based on degree centrality. The nodes denoted with ■ were only selected by the GA 
and not the degree strategy. Symbol □ denotes the nodes that were only selected by the degree strategy. 



4. Provide protection for the nodes found be the genetic algorithm, see section 2.5 



In order to keep the comparison fair, each strategy is allowed to protect exactly k nodes. The parameters 
k and / are set to 10 and 100, respectively, for the high school network and to 50 and 200 for the air 
transportation network (see Table [j}. 

After applying a strategy on a network, we run 500 SIR simulations and count the number of nodes 
that got infected by the virus. Before each SIR run we randomly select a single node to be infectious (the 
others are susceptible) in order to account for the fact that one is normally unaware where the virus starts. 
The results are presented in Fig. [3j The 'no protection' case represents the situation when no preparations 
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Fig. 5: The degree distributions of the high school and air transportation networks. The bars marked 
black denote the number of nodes with that particular degree that were targeted by the genetic algorithm. 
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were undertaken. For both networks this would have disastrous consequences. In the air transportation 
network, almost all nodes will definitely get infected. The genetic algorithm and the degree centrality 
approach seem to be preferred for both networks. The GA seems to outperform the degree strategy, but 
only slightly. The methods based on betweenness and eigenvector centralities do improve the situation, 
but are clearly not optimal. Note that the performance of these two approaches switch places between 
the graphs. The reason why eigenvector centralities might outperform the betweenness approach when 
applied to the air transportation network might lie in the fact that degree and eigenvector centrality 
expose a higher correlation in this graph than in the high school network, see Fig. [2}) and[2^. 

The genetic algorithm outperforms the strategy based on degree centrality only in the case of the 
US air transportation network. For the high school network their performances are comparable. Fig. |4] 
presents the topology of both networks and the differences between the sets of nodes proposed by the 
GA and the degree centrality strategy. The nodes depicted with A are the ones that were suggested by 
both strategies. The nodes depicted with ■ were only selected by the GA and not the degree strategy. 
Symbol □ denotes the nodes that were only selected by the degree strategy. Note that the GA manages 
to locate nodes that when protected will help to 'separate' clusters, i.e., it becomes harder for the virus 
to move from one cluster to the next. 

Fig. [5] depicts the degree distributions of both the high school and the air transportation network. 
The bars marked black denote the number of nodes with that particular degree that were targeted by the 
genetic algorithm. Note that the nodes selected by the GA are not all in the tail but also appear at the 
beginning and in the middle of the degree distribution. 

4 Conclusions & Discussion 

The work presented in this paper was undertaken to design and evaluate the effectiveness of a heuristic 
optimization method for virus spread inhibition, which, in contrast to earlier research ||7l[8l [T0l[TTl[T2l l9l. 
does not rely explicitly on any centrality measures during its search. 

We found that the genetic algorithm proposed in this paper seems to outperform the standard ap- 
proaches in the literature. This is a remarkable feat since the genetic algorithm searches through the 
immense set of possible subsets of nodes and does not base its decision directly on various centrality 
rankings. Fig. [4] shows that the genetic algorithm is able to find an adequate subset of nodes that need 
to be protected: it is able to construct a set of nodes that not only scores high on degree, but also on 
betweenness centrality. The virus is, thus, restricted in its movements within a cluster (degree) but also 
from one cluster to the other (betweenness). 

The degree distributions in Fig.[5]show the interesting result that the nodes with the biggest influence 
on the spread of the virus over the network are not necessarily to be found in the tail of the degree 
distribution; among them there are the nodes with no high degree at all, still, by being immunized, 
they seem to reduce the number of infected significantly. Future research is required to identify why 
exactly these nodes are of importance, and, most importantly, whether there are ways to simplify the 
identification of these individuals. Perhaps (non-topological) node attributes can provide important clues 
for the selection of individuals to be immunized? Being able to locate these individuals effectively 
without the need for the full topology of the network will aid tremendously in mitigating the impact of 
the virus on a population. 

The developed method could be applied rather easily to a wide variety of situations that require 
urgent computing since the method is not limited to a specific type of network: (near-)optimal solutions 
can be found for both regular and heterogenous networks. In addition, one does not have to restrict to 
one particular centrality measure. The search space of the genetic algorithm can be seamlessly extended 
for various fields of applications, regardless of the type of metrics used. 

The presented method is recommended for topologically heterogenous networks that are character- 
ized by the absence of strong correlations between the various node centralities and the presence of 
modularity. The authors expect that for homogenous networks the traditional degree centrality approach 
(i.e., immunizing the nodes with the highest degree) is to be preferred. In addition, the GA requires 
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more computation time. Especially when time is pressing, one is, therefore, recommended to resort to 
the traditional and fast-to-compute strategies. 

Further research might be done to explore different sets of metrics and heuristics in order to widen 
or narrow the search space. In addition, it might be interesting to take edge centralities into account 
as well, although the authors expect that this would not make a significant difference, since edge and 
node centralities are often highly correlated. The parameters for the genetic algorithm were set rather 
intuitively. Better results might be obtained when different parameter settings are applied. In this paper 
we only explored one particular compartmental model, the SIR model (see section [2~4] ). The method 
can be extended to different kind of models as well HI [3. It would be interesting to see how the method 
proposed here would perform when different virus dynamics are taking into account. The applicability 
and performance of other optimization methods such as simulated annealing might be a subject of future 
research as well. 

This paper has shown the promise of applying a genetic algorithm in order to mitigate the impact 
of a virus attack on a population. The authors would like to express their hope that continuing research 
along the lines set out in this paper will assist in gaining a better understanding of virus inhibition, and 
will ultimately provide better ways to avert future epidemics. 
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